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Abstract 

A measurement is made of the cross section for the process ep eXY in deep- 
inelastic scattering with the HI detector at HERA. The cross section is presented 
in terms of a differential structure function F^^^\xp, f3, Q"^) of the proton over the 
kinematic range 4.5 < < 75GeV^. The dependence of F^^^^ on Xp is found to 
vary with /3, demonstrating that a factorisation of F^^^^ with a single diffractive 
flux independent of /? and is not tenable. An interpretation in which a leading 
diffractive exchange and a subleading reggeon contribute to F^^^^ reproduces well 
the Xjp dependence of F^^^^ with values for the pomeron and subleading reggeon 
intercepts of 0^(0) = 1.203 ± 0.020(stat.) ± 0.013(sys.)t[{:|j^g (model) and 0^ (0) = 
0.50 ± 0.11(stat.) ± 0.11(sys.)taio(™odel), respectively. A fit is performed of the 
data using a QCD motivated model, in which parton distributions are assigned to 
the leading and subleading exchanges. In this model, the majority of the momentum 
of the pomeron must be carried by gluons in order for the data to be well described. 
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1 Introduction 



For a complete understanding of the strong interaction it is necessary not only to sur- 
vey short range, hard partonic interactions, but also the long range soft interactions that 
make up the bulk of hadron-hadron and photon-hadron cross sections at high energy. The 
lack of a hard scale in such soft processes makes calculations within quantum chromody- 
namic theory (QCD) problematic, and so these processes are usually described using the 
phenomenology of Regge theory [|l| by the t-channel exchange of mesons and, at high 
energy, by the leading vacuum singularity, the pomeronQ (P) [Q. Measurements of deep- 
inelastic diffraction, made possible for the first time with the HERA electron-proton (ep) 
collider^, offer a unique chance to study the soft scattering process with a hard virtual 
photon probe. It is hoped that precision measurements of diffraction in deep-inelastic 
scattering (DIS) will guide the formulation of a 'reggeon field theory' that encompasses 
both hard and soft aspects of the strong force |^,^ and will therefore provide a more 
complete understanding of QCD. 

Following the observation of deep-inelastic scattering events with a rapidity gap 
first measurements were made of the contribution of such events to the DIS cross section, 
quantified in terms of the differential structure function F^^^^ Hii]- These measurements 
demonstrated that the rapidity gap events could be attributed to an ep process in which 
the interaction is dominantly diffractive. No substantial dependence of F^^^^ was 
observed, thus indicating that the substructure of diffractive exchange was of a point-like, 
presumably partonic, nature. 

Presented here is a new measurement of F^^^\ based upon the data collected by 
HI during 1994. The measurement is made using data in which a forward rapidity gap 
is observed so that the hadronic final state is partitioned into two distinct subsystems 
X and Y, where Y is the system with rapidity closest to that of the incident proton. 
System X is measured with the main components of the HI detector, whereas system 
Y is not detected directly but the location of the rapidity gap implies its mass is less 
than 1.6 GeV. The data amount to a factor ~ 10 increase in statistics compared to our 
previous analysis [§] and they permit a substantially more accurate measurement of F.[^^^^ 
over a wider kinematic range. Further, by utilising data taken with the ep interaction 
vertex shifted by ~ 70 cm along the proton beam direction, the measurement is extended 
to considerably lower values of than hitherto possible. 

2 Cross Section Definition and Kinematic Variables 

In the analysis presented in this paper the cross section is defined in terms of a topological 
decomposition of the final state hadrons. The hadronic final state, defined as all final state 
particles excluding the scattered electron, is separated into two distinct systems labelled 
X and Y. These are taken to be the two hadronic systems which are separated by the 
largest interval in rapidity between final state hadrons, where rapidity is calculated in the 

^In this paper the terms 'pomeron exchange' and 'diffraction' are used synonymously. 
^HERA is able to collide both positrons and electrons with protons. In this paper the term 'electron' 
is used to describe generically electrons or positrons. 
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virtual photon-proton centre of mass frame. The system closer in rapidity to that of the 
incident proton is taken to be the subsystem Y. A schematic diagram of an ep scattering, 
with this decompostion of the final state, is shown in figure |l|. 

The kinematics of the process can be described by the invariant masses of the two 
subsystems, and M^, and the Lorentz scalars 

a:=^, y=^, Q' = -q', t = {P-pyf. (1) 
IF ■ q F ■ k 

Here P and k are the 4-momenta of the incident proton and electron respectively, py is the 
4-momentum of subsystem Y and q is the 4-momentum of the exchanged virtual photon 
(7*) coupling to the electron. The quantities x, y and Q"^ are the usual variables used 
to describe DIS; t is the square of the 4-momentum transferred at the proton-F vertex. 
Further useful variables are 

B= (2) 

^ 2q-{F-py) Q^ + Ml-f 

Here = {q + Ff is the square of the centre of mass energy of the virtual photon-proton 
system and Mp is the proton mass. 

The decomposition of the final state has been chosen to most clearly differentiate 
between various exchange mechanisms. Colourless exchanges can produce states in which 
X and Y are clearly separated in rapidity and both and My are small compared 
with W. For colour exchange processes, however, the location of the gap is defined by 
random fluctuations, the average gap size is small and one or both of and is 
generally large. In this paper events are selected (see section such that x^^ < 0.05 
and My < 1.6 GeV. This ensures that both <^ W and My < W, so that X and Y 
are well separated and the cross section is dominated by colourless exchange processes. In 
the infinite momentum frame of the incident proton, x^ may be interpreted as the ratio 
of the momentum carried by the colourless exchange to that of the incident proton, and 
/3 as the ratio of the momentum carried by the quark coupling to the virtual boson to 
that of the colourless exchange. 

In principle it is possible to measure a differential cross section with events decomposed 
in the scheme discussed above as a function of five independent variables, and to define a 
differential structure function F^^^\xp, jSjQ"^, M^jt) 



dxd/3dQ2dMy dt \^ ^ 2(l + i?«(5))^ 

Here is the fine structure constant and Q^, a;^, t. My) is the ratio of the 

longitudinal to the transverse photon cross sections. For the data used in this paper, 
however, the system Y is not measured directly, so that a differential measurement as 
a function of My and t is not possible. Instead a measurement of the triple differential 
cross section dld/Tdg^^ made, where an implicit integration is performed over the two 
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unmeasured quantities such that My < 1.6 GeV and |t| < 1 GeV^. The cross section is 
chosen to include, and is expected to be dominated by, the single dissociative process 



7*p — > Xp. The differential structure function F2^'^'^'*(Xp, /?, Q^) is defined and, following 
the usual convention is extracted with the assumtion of R^^^"* = 0. 



d^ 



dxd/?dQ2 (32Q4 



(5) 



The influence of a non-zero value of i?^*^^^ is discussed further in section ^. 



3 Simulation of Deep-inelastic ep Interactions 

In this analysis Monte Carlo simulations are used to correct the data for the effects of 
losses and smearing. A mixture of several models is used in order to best reproduce the 
data since no single event generator gives a complete description of all possible DIS final 
states. Each generated event sample includes a full simulation of the HI detector and is 
passed through the same analysis chain as the data. 

The RAPGAP generator ||T^ is used to generate the process ep —>■ eXp in the kine- 
matic region > 1.1 GeV and < 0.1. The model has two components in which the 
proton couples to a pomeron (P) or a meson (IR), which subsequently undergoes a hard 
collision with the incident electron. The pomeron and meson fluxes are taken to follow a 
simple Regge motivated form with an exponential t dependence so that 

B t B t 

F^^'\x^,(3,Q\t) = ^^4fT<(AQ') + -^^F^i/3,Q% (6) 

where the pomeron and meson trajectories are assumed to have the linear form, a{t) = 
a(0) + a't. The values for the intercepts ttjp(O) = 1.18 and «if^(0) = 0.6 are chosen such 
that F^^^\xjp, PjQ"^) from the generator is in close agreement with the results of the 
phenomenological fit presented in section ^ The values for the slope parameters, 
and Bj^, and the shrinkage parameters, a'^ and a'^, are not constrained by this analysis. 
The value assumed for the slope parameter of the pomeron, B^p = 6 GeV~^, matches 
the value obtained in a DIS measurement in which the leading proton was detected ||13|| . 
Since the shrinkage in the forward peak of the diffractive cross section has not been 
measured at high Q"^ and, since the value used in the generator leads to negligible changes 
to the measurements presented in this paper, a'^ is set to zero. The assumed values of 
B^ = 2 GeV-2 and a' = 0.9 GeV'^ are taken from hadron scattering experiments |H| , P^ 



The pomeron structure function Q^) assumed in the generator is derived from 

parton distributions taken from a QCD fit to the present data similar to that described 
in section |^. The meson structure function F2 {(3, Q"^) is taken from the GRV parameter- 
isation of the structure function of the pion |TB|. The structure functions are assumed to 
evolve with according to the leading order DGLAP evolution equations [|l^. In addi- 
tion to the 0{aem) quark-parton model diagram (eg — >• eg), boson-gluon fusion {eg — > egg) 
and QCD-Compton (eg — >• eqg) processes are generated according to the 0[aemOis) matrix 
elements. Higher order QCD corrections are provided either by the colour dipole model 
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as implemented in ARIADNE [|r^ or parton showers similar to those used in simulations 
of inclusive DIS (for example see LEPTO |]19|)- QED radiative processes are included via 
an interface to the program HERACLES |20[ . 

The low mass region of diffractive production {M^ < 1.1 GeV) is modelled with the 
DIFFVM ||2l|] generator, which simulates the production of the vector mesons p(770), 
07(782) and 0(1020). The cross sections are taken from measurements for p(770) and 
0(1020) production in DIS [|1-|2| and a ratio of c<;(782) : p(770) of 1 : 9 as given by 
quark counting rules. The and W dependences of the cross sections are assumed to 
follow those measured for p(770) production [p2| . 

The DIFFVM generator is also used to model proton dissociation. The cross section 
for proton dissociation is assumed to have a dcr/dM^ oc 1/M^ dependence at large M^, 



with a parameterisation based on measurements of proton dissociation on deuterium 
for the low mass region. Two different decay mechanisms are considered: decay with 
limited pt following the Lund fragmentation scheme ||2^, or isotropic phase space decay 



with multiplicity following KNO scaling laws The ratio of proton dissociation to 

elastic proton interactions is assumed to be 1. 

The kinematic region of high My or high (My > 5 GeV or > 0.1) is mod- 
elled with DJANGO [^. This model is a deep-inelastic ep generator based on the QCD 
substructure of the proton, with a colour string between the proton remnant and the 
current fragmentation region. Higher order QCD corrections are provided by the colour 
dipole model and QED radiative effects are implemented via an interface to HERACLES. 
DJANGO has been shown to reproduce well many aspects of the DIS final state in- 
cluding measurements of energy flows in the forward region 
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The normalisation for 

DJANGO was chosen such that the sum of the cross sections from RAPGAP, DIFFVM 
and DJANGO is in approximate agreement with that corresponding to the measured 
proton structure function F2{x,Q'^) 



The PHOJET |22[ generator was used to estimate background from photoproduction 
events. In PHOJET events are generated as five distinct subprocesses: non-diffraction, 
elastic vector meson production, vector meson production with proton dissociation, single 
photon diffraction dissociation and double diffraction dissociation. These subprocesses are 
mixed such that the total photoproduction cross section and the subprocess ratios are in 
agreement with measurements made at HERA [Q. 

Background events arising from the QED Compton process {ep epry) were simulated 



using the COMPOS generator [34 



4 Experimental Technique 

The data presented here were obtained by HI in 1994, when HERA was colliding 27.5 GeV 
positrons with protons at 820 GeV. The results shown are based on an integrated lumi- 
nosity of 1.96 pb^^ with the interaction point in its usual position relative to HI. A 
further 0.06 pb^^ of data, taken with the interaction point shifted by 70 cm in the proton 
direction, are also used. These data enable measurements to be made at lower values of 
Q^. The two event samples are referred to as the nominal vertex and shifted vertex data 
respectively. 
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4.1 HI Apparatus 



The HI detector is a composite, multi-purpose apparatus, designed to investigate the final 
states of high energy electron-proton interactions. Given here is only a brief outline of the 
components of the HI detector of relevance to the present analysis; for a more detailed 



technical description see I^Sj. A right handed co-ordinate system is employed at HI that 
has its z-axis pointing in the proton beam, or forward, direction. 

Surrounding the interaction point are the central and forward tracking detectors. 
These are situated in a uniform 1.15 T magnetic field, enabling momentum measurement 
of charged particles over the pseudorapidity range —1.5 < i] < 2.8.[| A multi-wire propor- 
tional chamber (BPC) gives additional acceptance for charged particles in the backward 
region (—3.0 < rj < —1.5). Surrounding the track detectors in the forward and cen- 
tral directions (— 1.5 < < 3.7) is a fine grained liquid argon calorimeter (LAr) and in 
the backward region (—3.5 < rj < —1.5) is a lead-scintillator electromagnetic calorime- 
ter (BEMC). The main calorimeters are surrounded by an instrumented iron backing 
calorimeter. 

Components of the HI detector positioned forward of the LAr are used in this analysis 
to tag particles at high values of pseudorapidity. These 'forward detectors' are sensitive to 
the secondaries produced by the rescattering of primary particles with the beam pipe and 
adjacent material, and are thus sensitive to larger pseudorapidities than would be expected 
on the basis of their geometrical acceptance alone. A copper-silicon plug calorimeter 
(PLUG) allows energy measurements to be made over the range 3.5 < rj < 5.5. Three 
double layers of drift chambers of the forward muon detector (FMD) are sensitive to 
particles produced at higher pseudorapidities 5.0 < ?7 < 6.5. The proton remnant tagger 
(PRT), consisting of seven double layers of scintillators, covers the region 6.0 < rj < 7.5. 

Two electromagnetic calorimeters (LUMI) situated downstream in the electron beam 
direction measure electrons and photons from the bremsstrahlung process ep epj for 
the purpose of luminosity determination. 



4.2 Event Selection 

Deep-inelastic scattering events are selected by requiring a compact energy deposit of at 
least 10.5 GeV in the BEMC and a BPC hit associated with this cluster. Background in 
the DIS sample is suppressed by requiring a vertex reconstructed by the track detectors 
within 30 cm of the mean z position of the vertex, a minimum and at least two recon- 
structed non-associated objects (tracks or clusters) besides the scattered electron in the 
final state. Further cuts, requiring an agreement between different kinematic reconstruc- 
tion methods, are introduced to reduce radiative corrections. The precise cuts employed 
are summarised in table 0a. A more extensive discussion of the DIS event selection at HI 



may be found in |31 



A subsample of events is selected where colourless exchange processes (such as diffrac- 
tion) are expected to dominate. This is achieved by requiring no activity above noise 
levels in any of the forward detectors or in the most forward part of the LAr (see table |I|b 

•^The pseudorapidity (?/ = — In tan 0/2) coverage of each detector component is given for the vertex in 
its nominal position. 
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(1) 


E'^ > 10.5 GeV 


(2) 


Q2 > 8.0 (3.5) GeV^ 


(3) 


156°(163°) < < 174°(176°) 


(4) 


El < 5 cm 


(5) 


£2 < 5 cm 


(6) 


kvtx - ^vtx < 30 cm 


(7) 


-^track + -^cluster > 1 


(8) 


-0.55 <y^-y^< 0.25 


(9) 


\y,,-yj< 0.2(0.3) 



b) 



(10) 


<ti < 3.3 


(11) 


^PLUG < 5 GeV 


(12) 


Nfmb < 2 


(13) 


A^PRT = 



Table 1 : Summary of event selection criteria for a) inclusive DIS events and b) additional cuts to select 
large rapidity gap events. The shifted vertex cuts are given in brackets where they are different from 
the nominal vertex cuts. Section 4.3 gives details of the methods used to reconstruct the DIS kinematic 
variables. Cuts 1-3 ensure the scattered electron is within the acceptance region of the BPC and BEMC; 
cuts 4 and 5 are electron identification estimators, ei is the electron cluster radius, £2 is the distance from 
the closest hit in the BPC to the centroid of the electron cluster ||3^]; cuts 6 and 7 require there to be 
at least 2 reconstructed objects in the hadronic final state of which one must be a track pointing to the 
interaction region; cuts 8 and 9 ensure agreement between different kinematic reconstruction techniques; 
cut 10 requires that there is no activity observed in LAr above noise levels for 77 > 3.3; cuts 11-13 ensure 
that there is no activity above noise levels in the forward detectors p&. 



for details of the cuts). The selection implies that systems X and Y (as defined in sec- 
tion]^) are well separated, such that X is fully contained in the main detector components 
and Y passes unobserved into the forward beam pipe. The maximum allowed value of 
pseudorapidity for particles from the system X of 3.3 (cut (10) in table |l|b) ensures that 
Xjp ^ 0.05. The fact that no part of the system Y is detected in the forward detectors 
restricts < 1.6 GeV and 1 1 1^ 1 GeV^. This sample of events is referred to hereafter 
as the 'large rapidity gap' event sample. The number of selected events is 13956 in the 
nominal vertex sample and 760 in the shifted vertex sample. 

The only backgrounds in the large rapidity gap event sample that were found to 
be significant were from photoproduction and the QED Compton process ep — > epj. 
Photoproduction background arises when the scattered electron escapes undetected along 
the beam pipe and a final state hadron fakes a BEMC electron cluster. This background 
was estimated using the PHOJET simulation and is found to contribute most significantly 
at high y and low Q^. The overall photoproduction background was estimated to be 0.3 % 
of the large rapidity gap event sample and at most 13 % for any single data bin. 

QED Compton events were studied with the COMPOS generator and found to be 
effectively suppressed by cut (7) in table |l|. Residual QED Compton background arises 
when the final state photon converts in material in front of the BEMC, producing a 
BPC hit and is mistaken for the scattered electron. The true electron and some residual 
calorimeter noise fake the hadronic final state. This background was found to contribute 
most significantly at high (3 and high y. The overall QED Compton background is esti- 
mated to be 0.6% of the large rapidity gap event sample and contributes a maximum of 
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15 % in any single bin. 

Each background source was subtracted bin by bin. 



4.3 Reconstruction of the Kinematic Variables 

The kinematic variables, defined in section ^, are calculated from the energy and the polar 
angle of the scattered electron {E'^, 6^) and from the 4-vector of the hadronic final state 
{Eh, Pxhi Pyhi Pzh)- E'^ is rccoustructed from the cluster with the largest energy in the 
BEMC passing the cuts described in table 0, and 6e is the angle defined by the z-vertex 
and the BPC hit associated with the electron. The hadronic 4-vector is reconstructed 
by combining information from the track detectors and calorimeters using an energy fiow 
algorithm, which improves the resolution compared to methods relying on single detector 
components [p7|] . 

A key feature of the HI experiment is its ability to measure both the electron and 
hadronic final state well, thus providing a number of methods for reconstructing the 
kinematic variables. In this analysis the kinematics are reconstructed using a method 
that combines the electron method, in which the kinematics are reconstructed solely 
from the scattered electron, and the double angle method, which uses only angles of the 



hadronic final state and the scattered electron The two methods are combined in 
order to utilise the optimal resolution of the electron method at high y and the double 
angle method at low y. We reconstruct y as 

y = yl + yDAiX-yDA)i (7) 

where y =l-^.u.\ej2) and y,^ = ^^"^(l + cosg) ^ 

^ ' ' sin7 + sin^e + sin(0e + 7) 

•^v, Plh+Plh-{Eh-p,h? 
with cos7= J (9) 

Pxh + Pyh + [Eh - PzhV 

and Ef, is the electron beam energy. As demonstrated by figure 0(a), this '?/ weighted 
averaging method' provides excellent resolution in y over its full range. is reconstructed 
from y and the angle of the scattered electron 

^ ^EKi-y) .,0) 

^ tan2(e,/2) • 

Two additional kinematic reconstruction methods are used to reconstruct y and/or 
for the purposes of cross checks and event selection cuts. These are the hadron method 



for which y,^ = {Eh — pzh)l{'2Ef,) and the S method pO| for which y^ = y^J {l-\-y^^—y^) and 
Q'^ = 4:E'^ siia^^Oe)/ {1 — y^)- The Bjorken scaling variable x is calculated in all methods 
as Q'^/sy where s is the ep centre of mass energy. 

For the large rapidity gap event sample the system X is fully contained within the 
main detector, so that the visible hadronic final state corresponds to the system X. It is 
thus possible to reconstruct the variables , f5 and well. The invariant mass is 
reconstructed with a method that utilises the optimised reconstruction of y 

mI = {eI-pI,-pI,-pI,)^. (11) 
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A constant scale factor of 1.10, determined from the simulations, is applied to the re- 
constructed to correct for residual losses, mainly due to dead material in front of 
the calorimeters. As demonstrated in figure |(b), the resulting is well reconstructed 
across the kinematic range of the measurement with a resolution of about 25 %. The other 
kinematic quantities are determined using equations (0) and (|]), with the approximation 

The accuracy of the reconstruction method and the ability of the simulation to describe 
the data were checked in a variety of ways. Figure ^(a) shows the ratio of the transverse 
momentum of the visible hadrons to that of the electron, which for a perfect detector 
should have a narrow width and peak at 1. At the detector level it can be seen that 
the peak is smeared with a ~ 30%, but with a mean close to 1, showing that there 
are no substantial losses of final state hadrons. Figure ||(b) shows that there is good 
agreement between y when reconstructed using the electron or double angle methods 
before combination. For both these distributions the simulation is in good agreement with 
the data, showing that the resolutions and biases of the detector are well understood. 

Since an understanding of bin migrations depends not only on the resolution in the 
kinematic quantities, but also on the form of the differential cross section, it is important 
that the cross section used as the input to the simulations closely matches that of the 
data. Figures ||(c) and ||(d), which show the reconstructed [3 and 6e spectra, give an 
indication that this is indeed the case. Similar agreement is found for the x^, and y 
spectra. 

4.4 Evaluation of the Differential Structure Function 

The large rapidity gap event sample is used to measure the differential cross section 
dxd/3dQ2 for the process ep — > eXY as defined in section ^ for 4.5 < < 75 GeV^, 
0.0002 < < 0.04 and 0.04 < (3 < 0.9. The data are divided into three-dimensional 
bins in x, (3 and Q^. The Monte Carlo models are used to correct bin by bin for losses, 
smearing and backgrounds. No difference was observed within the statistical errors of 
the simulations between the acceptances when evaluated using the ARIADNE or parton 
showers fragmentation options within RAPGAP. The average of the two acceptances was 
therefore taken. An additional overall correction of 1.11 ± 0.05 is applied to account for 
events lost due to noise in the forward detectors, which is not simulated in the Monte 
Carlo models. This correction is evaluated from randomly triggered events, in which there 
is no requirement of particle activity in the detector. Radiative corrections are applied 
bin by bin using the HERACLES interface to RAPGAP. 

Since the cross section is expected to vary smoothly with the kinematic variables for 
bins that lie away from the resonance region of the light vector mesons (p'^, uj, (p, p'), the 
cross section is corrected to that at the bin centre (xc, Pc, Ql) for bins where Mxc > 2 GeV. 
For these bins the main deviation from a smooth variation of the cross section is expected 
to arise from elastic J/tp and ip{2S) production, which contribute a maximum to any 
single bin of approximately 5% and 1% respectively |^,^. These contributions are not 
subtracted but a systematic error equal to the expected magnitude of the contribution is 
added to the systematic error. For the bins that have a mass interval such that they cover 
one or more of the light vector mesons {Mxc < 2 GeV), the dependence of the cross section 
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on X and is expected to be smooth, but there are hkely to be substantial deviations 
from a smooth dependence of the cross section as a function of f3. For this reason the 
cross section is corrected to the bin centre only in (xc, Q^) and in f3 it corresponds to the 
average over the range covered by the bin. 

For all selected bins the acceptance is required to be above 0.4, and the purity (i.e. 
the fraction of the number of events both generated and reconstructed in a bin to the 
number reconstructed in the bin) is required to be greater than 0.2. The chosen bins have 
an average acceptance of 0.82 and an average purity of 0.45. The two data samples are 
combined by taking the shifted vertex data where there is no nominal vertex coverage 
{Q^ < 9 GeV^) and the nominal vertex data elsewhere {Q^ > 9 GeV^). 

The statistical errors are formed by adding in quadrature the contribution from the 
data and that due to the finite statistics of the Monte Carlo simulations. Systematic 
errors arise from the following sources: 

• The error due to uncertainties in the reconstruction of the scattered electron, es- 
timated by changing the electromagnetic energy scale of the BEMC by ±1% and 
shifting the electron polar angle by ±1 mrad. 

• Uncertainties in the reconstruction of the hadronic final state, estimated by changing 
the hadronic energy scale of LAr by ±5%, that of the BEMC by ±20% and by 
shifting the energy fraction carried by the tracks by ±3%. 

• Errors due to our a priori ignorance of the input structure function as used in 
the simulations, estimated by reweighting the generated events by factors 

and (1 ± 0.3/3). These factors are chosen such that resulting generated structure 
function is changed by an amount larger than the final errors on the data points. 
An additional error is assigned to account for smearing of events from kinematic 
regions where the differential structure function is not measured. This error is 
conservatively estimated by changing the number of events generated with > 0.1 
and/or M^. > 5 GeV by ±50%. 

• Errors due to the uncertainty in the t dependence, estimated by reweighting by e^^*. 

• The uncertainty in the correction due to smearing of events across the boundary 
My = 1.6 GeV, estimated by changing the PLUG energy scale by ±30%, by chang- 
ing the PRT efficiency by ±20%, by taking the difference between the correction 
evaluated using isotropic and limited pt proton dissociation models and by varying 
the ratio of double dissociation to single photon dissociation assumed in the sim- 
ulation between 0.5 and 2.0, such that the range covers that measured in various 
proton dissociation processes in DIS and photoproduction p5|,^,|42||. 



The uncertainty in the number of background events, estimated as 50% of the 
photoproduction and 30% of the QED Compton background. 

The uncertainty in the bin centre correction due to contributions from high mass 
resonances, estimated by assuming an error of 5% for all bins which contain the J/ip 
and 1% for all other bins to account for the ip{2S) and other possible contributions. 
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• An error of 2% per point is assumed to account for any detector inefficiencies not 
included in the simulations 

• An overall normalisation uncertainty of 4.8% for the nominal vertex data and 6.0% 
for shifted vertex data, evaluated by adding in quadrature a 4.5% error on the 
uncertainty due to the correction for events lost due to forward detector noise and 
the errors on the integrated luminosity of 1.5% and 3.9% for nominal and shifted 
vertex data. 

To check the stability of the results, the analysis was repeated by widely varying 
the experimental cuts and by changing kinematic reconstruction techniques. Figure ^(a) 
shows that there is reasonable agreement between the results obtained from the nominal 
and shifted vertex data samples for the standard analysis in the region of overlap. In 
figure ^(b) it can be seen that the results do not change significantly if the S method is 
used to reconstruct the kinematics. Figures |(c), (d) and (e) demonstrate that changing 
the rapidity gap requirement does not greatly change the results. In figures ||(f) it can 
be seen that even without the requirement of a rapidity gap, when the measurement is 
made solely on the basis of the reconstructed values in the inclusive data sample, the 
results obtained are well within the quoted systematic uncertainty. 

The measured values of F2^^^''(Xjp, /5, Q^), evaluated using equation (^, are listed in 
table 1^ and are also available from . 

The quantity x^F2^'^\x-p^ (i^Q"^) is shownQ in figure ^ as a function of for dif- 
ferent (3 and values. It can be seen that over most of the kinematic region covered 
a;pF2^*-^^(Xp, /3, Q^) is either falling, or approximately constant, as function of increasing 
Xjp. Qualitatively this is consistent with our earlier observation of a dominant diffrac- 
tive contribution with an intercept ^^(O) 1. However, the more precise data of the 
present analysis show that there is a clear change in the x^ dependence at low (3 where 
there is a tendency for it to increase at the highest x^. A quantitative study of this x^ 
behaviour is given in terms of Regge amplitudes in the following section. 



5 Regge Parameterisation of the Data 

An appropriate phenomenological framework to study the behaviour of the DIS colour 
singlet exchange cross section is provided by Regge theory |[l|,3. This theory has proved 
to be remarkably successful in describing high energy hadron-hadron and hadron-photon 
cross sections (see for example 0,^^^) as well as early measurements of F2 
r3| . Regge theory in its simple formQ, when applied to deep-inelastic scattering, provides 
information on the x^ dependence at fixed /3, and t. Since there is much theoretical 
uncertainty concerning the precise form of the Regge model which is appropriate for colour 
singlet exchange in DIS, several fits to the data are performed with different assumptions. 
In all fits the values are quoted for statistical errors only. 



'^x-pF2^^^ {x-p , (3, Q^) is plotted, rather than F^*-^-* (x^ , /3, Q^), in order to show more clearly any devi- 
ations of the data from a simple power law behaviour on x-^ . 

^No attempt has been made in this paper to describe F^^^\xjp, I3,Q'^) using 'triple-Regge' theory, 
which would provide an additional prediction for the dependence of the cross section on (3. 
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For fit A it is assumed that F-^'^^^ may be factored into a pomeron flux and pomeron 



structure function 111 



The pomeron flux is taken to follow a Regge behaviour with a peripheral t dependence 
and a linear trajectory a^pit) = ap{0) + a'^t such that 

where \tmin\ is the minimum kinematically allowed value of |t| and tcut = — 1 GeV^ is the 
limit of the measurement. The values for and a' cannot be constrained in the fit and 
are assumed to be the same as those obtained from analyses of hadron-hadron data. These 
values give consistent results with a measurement of the t dependence of the colour singlet 
exchange cross section averaged over the other kinematic quantities |jl3|. Table |^ lists the 



assumed values for all fixed parameters used in this and subsequent fits. The value of 
0^(0) is left as a free parameter in the fit as is (/?, Q^) at each where there 

are data. A complication arises from the lack of knowledge of R^^^\P, Q^, Xp,t, M^), the 
ratio of the longitudinal to the transverse cross section. The effect of a non-zero value for 
/j^(5) produces a decrease in the F^^^^^ extracted at large y relative to that obtained with 
the assumption R^^^^ = 0. The maximum value of y for any of the HI data is at ?/ = 0.55 
which leads to a maximum shift in F^^^^ of 34% if i?^*^^^ = oo. Since these shifts are 
large compared with the errors on the data points only those data with y < 0.45 are fitted 
thus limiting the shift to a maximum of 18%. The difference when R^^^^ = oo instead of 
Jl^i^) = is included as part of the 'model dependence' error for all fits. 



Quantity 


Value 


Source 




" -0 




< 


0.26 ±0.26 GeV-^ 




46 




a' 

m 


0.90 ±0.10 GeV-^ 
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4.6 1^:1 GeV-^ 




46 






2.0 ±2.0 GeV-^ 
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Table 2: Values assumed for the fixed parameters in the Regge fits and the sources from which they 
are taken. Each quantity is defined in the text. Contributions to the model dependence errors are formed 
by repeating the fits after separately varying each parameter by the quoted uncertainties. 

Fit A yields a rather poor x^/ndf of 257.8/168. Changing the form of the depen- 
dence or the values of the fixed parameters of the fit, whilst still retaining the factorisation 
hypothesis, cannot explain the data simultaneously at both high and low f3. The exchange 
of a single factorisable Regge trajectory in the i-channel does not, therefore, provide an 
acceptable description of the dependence of F^^^'' on x^. 

A natural explanation of the observed factorisation breaking is the presence of a sub- 
leading exchange in addition to the leading pomeron. Such an exchange is necessary to 
describe the energy dependence of total, elastic and dissociative photoproduction and 
hadron-hadron cross sections [14, 37,47 
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In fits B and C we assume a subleading exchange, made up from tlie set of reggeons 
whicli lie on the approximately degenerate trajectory a]f^(t) — 0.55 + 0.9t and which 
carry the quantum numbers of the p, u, a or f meson. It is assumed that each of these 
exchanges can be expressed as the product of a flux and a meson structure function akin to 
equation (|l^) and that the structure function of each exchange is the same. The situation 
is complicated by the possibility of interference between the terms which describe the 
pomeron and the / meson trajectory^ so that F^^^'^ takes the form 



+21 f (x, ) {(3, Q^)F^ {(3,Q^), ( 14) 



with 



The interference flux is derived from the pomeron and meson fluxes with the phase com- 
pletely specified by the signature factors of the pomeron and the / [^] in the form of 

ftmin fir \ e(-^ip+-^iR)*/2 

f M = l^^^ COS (^-K(t) - a^it)] j ^,^(,)+,^(,),, dt. (16) 

The quantity / denotes the degree of coherence and is assumed to take any value from 
0, when no / contributes to the subleading exchange, and up to 1 for / dominance of 
the subleading reggeon In these fits / is assumed independent of /? and Q^. The 

quantities 

F2 (/?,g2) andFa (/3,Q2) 

are left as free fit parameters and a constraint 



that F2 (/?, Q ) > is imposed. 

For fit B no interference is assumed by setting 1 = 0. The fit gives an acceptable 
X^/ndf of 120.7/121 and describes the deviations from factorisation of F-^'^^^ well. Values 
of ajp(O) = 1.200±0.017(stat.)±0.011(sys.)l[J;[]^^(model) and a^(0) = 0.57±0.14(stat.) ± 
0.11(sys.)to;o7(model) are obtained. The systematic errors on both extracted quantities 
are estimated by adding in quadrature the differences obtained from the central fit value 



after imposing each of the variations described in section |4.4| . The model dependence 
error is estimated in a similar manner by varying the values of the fixed parameters of 
the fit over the ranges listed in table |^. These ranges represent our best estimate of the 
likely variation in the present data after examining the spread in the values obtained in 
other scattering processes. 

In fit C maximal interference is assumed between the pomeron and the meson tra- 
jectory (/ = 1). A very similar x^/ndf of 120.4/121 is obtained, with a similar value 
of ttjp jo) = 1.206 ± 0.022(stat.) ± 0.013(sys.)t[J:[]i^g(model) and a lower value of ^^^(0) = 
0.44 ± 0.08(stat.) ± 0.07(sys.)l[J:|]^(model). The resuhs of this fit are shown in figure ||. 
It can be seen that at constant the meson component and the interference part fall 
much more steeply as (3 increases than does the pomeron component. 

Since fits B and C both give a good description of the data and have similar x^/ndf 
it is not possible to determine, with the present data, whether or not interference plays 

^Interference between the pomeron and any of the other subleading mesons trajectories is forbidden 
by C- and G-parity conservation. 
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a significant role in DIS colour singlet exchange cross sections. Further precision mea- 
surements at high values of are needed to clarify this matter. The results of the two 
fits do, however, yield compatible results for the intercepts of the pomeron and meson 
exchanges, and, since it is likely that the truth lies somewhere between the solutions with 
/ = and / = 1, a best estimate for the values of the intercepts is obtained by averaging 
the two results and including the difference between the two fits in the model dependent 
error. This procedure yields 

a^{0) = 1.203 ± 0.020 (stat.) ± 0.013 (sys.)+°:°^° (model) 
a^(0) = 0.50 ±0.11 (stat.) ±0.11 (sys.)+no (model). 

The value of 0^(0) found in this analysis is significantly larger than that obtained from 
analyses of soft hadronic cross sections [0, J8|. It also exceeds the value «ip(0) = 1.068 ± 



0.016(stat.) ± 0.022(sys.) ± 0.041(model), which was obtained from an HI analysis of the 
photon diffractive dissociation cross section at = |3^. The value of ^^^(O) agrees 



with the value of ~ 0.55 obtained in an analysis of total hadronic cross sections 

An attempt is made to see if there are any further deviations from factorisation of the 
pomeron flux within the two-reggeon model described by equation (0). Two separate 
fits are performed that follow B, but now 0:^(0) is allowed to vary as a function of either 
/3 or so that a^{0) = ai ± a2(/3 - 0.5) or a^{0) = ± a'g logio(Q2/10GeV^). Both give 
acceptable x^/ndf of 120.5/120 and 120.7/120 respectively, but within the errors both fits 
show that the data are compatible with no change in 0^(0) with (3 or Q^: a2 = 0.03 ± 
0.10(stat.) ±0.04(sys.)t[|:[j^(model) , a'^ = 0.004 ± 0.071(stat.) ± 0.014(sys.)t[j:E][j|(model). 
Similar results are obtained for fits performed with full interference. 

We conclude that the data may be described using a Regge model with two components 
in which a pomeron contributes along with a subleading exchange. In this two reggeon 
scenario there is, within experimental errors, no evidence for a change in the pomeron 
intercept with (3 or over the kinematic region of this measurement. 



6 Deep-inelastic Structure of Diffractive Exchange 

6.1 QCD Extension of the Phenomenological Model 

In the universal factorisation hypothesis discussed in section ^ the structure functions 
F^{P, Q^) and F^{(3, Q^) describe the deep-inelastic structure of the pomeron and meson 
exchanges respectively. It has been suggested that the evolution of these structure 
functions may be understood in terms of parton dynamics and therefore perturbative 



QCD [11,49,5C 



To establish whether the data are consistent with such a hypothesis, the phenomeno- 
logical model presented in section | is further developed by assigning parton distribution 
functions to the pomeron and to the meson. Currently, there is no clear theoretical con- 
sensus on how colour singlet exchange processes may be understood within the context of 



QCD, perturbative or otherwise Therefore, a simple prescription is adopted in which 
the parton distributions of both the pomeron and the meson are parameterised in terms 
of non-perturbative input distributions at some low scale Qg. The parton distributions 
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of the pomeron and of the meson are evolved separately with increasing according to 



the DGLAP evolution equations. 

For the pomeron, a quark flavour singlet distribution {zj^q{z, Q"^) = u+u+d+d+s+s) 
and a gluon distribution {zj^g{z, Q"^)) are parameterised in terms of the coefficients Cj at 
a starting scale of = 3 GeV^ such that 



zTi{z, Q'^ = Q 



Y^C]-P,(2z-V, 



2 



■e^ (17) 



where z = Xijp is the fractional momentum of the pomeron carried by the struck parton. 
If the photon couples directly to a quark {i = q) then z = (3, whilst if the photon 
interacts with a gluon [i = g) then z is the fractional momentum carried by the gluon 
and < (3 < z. Pj{C) is the j*^ member in a set of Chebychev polynomials, which 
are chosen such that Pi = 1, P2 = C and Pj+i(C) = 2(Pj{Q — Pj^i{Q. A sum of n 
orthonormal polynomials is used so that the input distributions are free to adopt the 
widest possible range of forms for a given number of parameters. Any bias towards a 
particular solution due to the choice of the functional form of the input distribution is 
therefore minimised, and the number of terms included in the parameterisation may be 
matched to the significance of the data. The series is squared to ensure a positive definite 
parameterisation for all values of z and C*. 

For the evolution equations to be soluble, the parton distribution functions must tend 
to zero as z — >■ 1. This is achieved by introducing the exponential term with a positive 
value of the parameter a. Unless otherwise indicated, in the following fits a is set to 0.01 
such that this term only infiuences the parameterisation in the region z > 0.9. 

The functions zJ^i{z,Q'^) are evolved to higher using the next-to-leading order 
DGLAP evolution equations and the contribution to Ff{(5, Q"^) from charm quarks is 
calculated in the fixed fiavour scheme using the photon-gluon fusion prescription given 
in The contribution from heavier quarks is neglected. 

In the Regge model phenomenological fits B and C described in section ^ the fitted 
values of P2 (/5,g2) 

increase with decreasing [3 such that they are consistent in shape with 
parameterisations of the pion structure function. Since the errors on the current data 
are insensitive to any evolution of F2 Q^), the structure function of the subleading 
exchange is taken from a parameterisation of the pion multiplied by a single coefficient 
CjR to be determined from the data. 

No momentum sum rules are imposed because of the theoretical uncertainty in speci- 
fying the normalisation of the pomeron or meson fiuxes (/^/^(a^jp), and because 
it is not clear that such a sum rule is appropriate for the parton distributions of a virtual 
exchange. 



6.2 QCD Fits of F2^*'^(cc^,/3,Q2) 

In the following fits the values of F2^'^\l3^Q'^-,X-^) are calculated from F^{(3,Q'^) and 
F^{[3, Q"^) using equation (0) assuming no interference between the leading and sub- 
leading trajectories. All QCD fits are performed using the total statistical errors. As 
in the phenomenological fits, the data for y > 0.45 are excluded from all of the QCD 
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fits to limit the uncertainty introduced by the lack of knowledge of R^^^\ Data points 
corresponding to values of Mx less than 2 GeV are also excluded because of the large 
contribution to the cross section from the resonant production of vector mesons and other 
possible higher twist contributions in this region. 

For fit 1 only quarks are assumed to contribute to the structure of the pomeron at 
the starting scale of Ql = 3GeV^, and the values of 0^(0) and ^^^(O) are fixed to those 
determined in the Regge fit without interference (fit B). The parameters to be fitted are 
then Cj and the normalisation of the subleading component Cjr. Introducing successive 
terms from the Chebychev series, a x^/ndf of 314.2/159 is reached for 3 terms. If an 
additional term is introduced then the corresponding coefficient is consistent with 0, with 
no improvement of the quality of the fit. The resulting parameterisation of ^2^*^^'' is shown 
at a constant value of = 0.003 as a function of for different values of (3 in figure 
11(a), and as a function of (3 at fixed in figure 0(a). Also shown are the data for x^F.^^^^^ 
interpolated to x^ = 0.003 using fit B. The large x^/ndf for fit 1 arises from the failure 
to reproduce the scaling violations of F^^^\x-p, (3, Q^), which rises with increasing for 
(3 < 0.9. Therefore, the data cannot be described by a parameterisation in which the 
pomeron contains only quarks at = 3 GeV^. 

In fit 2 the first term in the expansion of the gluon distribution is introduced as 
an additional fit parameter, and a significantly better description of the data is obtained 
(x^/ndf = 186.7/156). The description of the data improves further with the introduction 
of the next two terms in the series to give x^/ndf = 176.3/154 for n = 3 (fit 3). This fit is 
shown in figures ^b) and|^(b) where it can be seen that the dependence of F^^^\xp, (3, Q'^) 
on and /5 is well reproduced. Introducing additional terms for the gluon distribution 
does not result in any further improvement of the quality of the fit. 

The parton distributions resulting from fits 2 and 3 are shown in figure |^ where both 
the sum of the light quark distributions and the gluon distribution are shown for each fit 
at the lowest and highest values included in the analysis as well as an intermediate 
value of = 12 GeV^. In both fits 2 and 3 the fraction of the momentum of the pomeron 
carried by gluons decreases with increasing from ~ 90% at = 4.5 GeV^ to ~ 80% 
at = 75 GeV^. A large gluon component in the pomeron is, therefore, favoured by the 
data. 

Repeating fit 3 with 0^(0) and «j^(0) as free parameters results in no change in either 
the coefficients Cf^ or the quality of the fit, and the values for 0^(0) and a^^^O) are 
unchanged from those obtained in the Regge fit B. The parameters a and Qq may be 
varied in fit 3 within the ranges < a < 0.1 and 2 < Qq < 4GeV^ without producing 
a change of more than one unit in the x^y aiid if a is allowed to vary in the fit then the 
value obtained is consistent with 0.01. A final fit is performed to assess the influence of 
the particular choice of the meson structure function by reweighting the parameterisation 
of by a factor (1 — /?)'''«. The value of 7jr obtained is consistent with zero. 

6.3 Discussion of the QCD Fits 

We conclude that it is possible to describe the data for F^^^^ in the region Mx > 2 GeV 
with a factorisable model in which the pomeron structure function evolves with ac- 
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cording to the NLO DGLAP evolution equations. Within the context of this model a 
large gluonic content of the pomeron is required, with a hard distribution in fractional 
momentum in order to achieve an adequate description of the data. The data exhibit 
some preference for a gluon distribution that is large at high z such that for low values of 
the majority of the momentum of the pomeron is carried by a single gluon, although 
the precise shape of the gluon distribution is not well constrained. It should be noted, 
however, that there are some unresolved theoretical problems associated with the use of 
the DGLAP evolution equations at very large z due to resummation effects |^^, which 
may be particularly important when the parton densities near z = 1 are very large as in 
fit 3. 

The fraction of the momentum of the pomeron carried by gluons is similar in both 
fits 2 and 3, from ~ 90% at = 4.5 GeV^ to ~ 80% at = 75 GeV^. Thus the data 
lend support to pictures of diffraction in DIS in which the dominant mechanism is boson- 
gluon fusion with the incoming gluon carrying a large fraction of the momentum of the 
pomeron 

Measurements of the properties of hadronic final states produced in diffractive hard 
interactions are essential to clarify the extent to which this partonic interpretation may 
be applied universally. For example, a large gluon distribution in the pomeron, peaked 
near z = 1, was found to describe well the measurements of diffractive dijet production 



in pp interactions [Q. The large gluon distribution favoured by the present analysis 
predicts a sizable cross section for diffractive dijet and charm production in both pp 
and ep interactions. To facilitate comparisons between such measurements [^] and the 
parton distribution functions presented here, parameterisations of the parton distributions 
functions for the NLO fits 1, 2 and 3 are available |H3|. Parameterisations are also available 



of the results of three leading order QCD fits which are identical to the NLO fits 1, 2 and 
3 in all other respects. The results of these additional fits, not shown here, lead to the 
same conclusions as those outlined above. 

The QCD analysis of F^^^^^ presented in this paper does not attempt to address the 
physical origin within QCD of either the intercepts of the leading and subleading ex- 
changes, or of the parton distributions at the starting scale of the perturbative evolution. 
It is hoped that a more complete understanding of colour singlet exchange interactions 
within QCD will establish the extent to which the DGLAP evolution prescription is ap- 
propriate, and also provide predictions for 0:^(0) and a^^iO). 



7 Summary 

The differential structure function F^^^\xp, PjQ"^) has been measured using the data 
taken in 1994 with the HI detector with greater accuracy and over a greater kinematic 
range than previously. At this level of precision it is found that ^2^*^^^ is inconsistent with 
a universal dependence on x^. A natural explanation of this breaking of factorisation is 
presented in which a leading pomeron exchange is accompanied by a subleading meson 
exchange. Fits to the data in which such a subleading term is considered yield an inter- 
cept for the pomeron of a^{0) = 1.203 ± 0.020(stat.) ± 0.013(sys.)l[]:J]i^g (model), which is 
significantly higher than the value obtained in similar parameterisations of soft hadronic 
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cross sections. No evidence is found for any variation of a^p (0) with or P in the range 
4.5 < (5^ < 75 GeV^ and 0.04 < P < 0.9. The value obtained for the intercept for the 
subleading exchange, a^{0) = 0.50 ± 0.11(stat.) ± 0.11(sys.)lg:?^(model), was found to 
agree well with values obtained from soft hadronic scattering processes. 

A model in which the deep-inelastic structure of the leading exchange is described by 
parton distributions which evolve according to the NLO DGLAP evolution equations is 
consistent with the data for F^^^\xjp, PjQ"^) in the region Mx > 2 GeV. Under such 
a hypothesis, the data require approximately 90% and 80% of the momentum of the 
pomeron to be carried by gluons at — 4.5 GeV^ and = 75 GeV^ respectively. A 
parametrisation which assumes only quarks in the diffractive exchange at — 3 GeV^ is 
excluded. 
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Figure 1: Schematic of an DIS scattering event ep — > eXY. 
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Figure 2: The correlation between the generated and reconstructed values of (a) y (reconstructed using 
the weighted average method) and (b) for simulated Monte Carlo mixture at nominal vertex. 
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Figure 3: The solid points show (a) the ratio of the transverse momentum of the hadrons to that of 
the electron, (b) the difference between y calculated from the electron and double angle methods, (c) the 

(3 distribution and (d) the electron polar angle distribution. All distributions are shown for the rapidity 
gap DIS event sample. Plots (a), (b) and (c) arc for the nominal vortex data and plot (d) is for the 
shifted vertex data. None of the distributions have been corrected for detector effects. The predictions 
from the simulation are shown for each plot as a solid histogram. 
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Figure 4: Effect of changing kinematic reconstruction techniques and rapidity gap selection criteria 
on the measured values of x^F^*-^-* at an example value of 18GeV^. The solid circles show the data 
measured using the standard analysis. In each of the six columns, these data are compared to those 
obtained with a modified analysis technique in which, (a) the shifted vertex data are compared (b) x 
and were reconstructed with the E method and solely from the hadronic final state, (c) cut 10 
(see table ||) is changed to TyJ^^x < 3, (d) cut 10 is not used, (e) cut 10 is replaced by rj^^l, < 2 and 
there is no requirement on the forward detectors (i.e. cuts 11-13 are not applied), (f) there is no rapidity 
gap selection at all (i.e. cuts 10-13 are not applied). The comparisons are plotted only for those bins 
satisfying the minimum acceptance and purity limits as for the standard analysis. 
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Figure 5: The solid points show the measured differential structure function plotted as 
XjpF^^^^ (Xjp , /3, Q^) against x^p for various f3 and values. The inner error bars represent the statistical 
errors and the outer error bars represent the statistical and systematic errors added in quadrature. The 
curves show the results of the Regge fit with interference (fit C) described in the text. The dashed curve 
is the contribution from the pomeron alone, the dotted curve is the pomeron plus interference and the 
continuous curve is the total contribution. The data that lie at values of x^p smaller than the extent of 
the curves have y > 0.45 and are excluded from the fit. 
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Figure 6: The quantity a;pJ'2^*^\a;jp,/3,(5^) at x^^ = 0.003 as a function of for different values of /3. 
In both (a) and (b) the solid points show the same data, interpolated to x^p = 0.003 using fit B (see text). 
In (a) the result of QCD fit 1 (in which only quarks contribute to the pomcron structure at = 3 GeV^) 
is superimposed. In (b) the preferred QCD fit 3 (in which both quarks and gluons contribute) is shown. 
The results of the QCD fits are only shown in the kinematic range used in the fit (corresponding to data 
for which Mx > 2GeV). In both figures the solid line represents the value of F^^^\xjp, P,Q'^), whilst 
the dotted line shows the contribution from the pomeron term only. 
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(b) Quarks + Gluons (Fit 3) 
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Figure 7: The quantity a;pi^2 (Xf^P^Q^) sA, = 0.003 as a function of /3 for different values of Q^. 
In both (a) and (b) the solid points show the same data, interpolated to x^p = 0.003 using fit B (see text). 
In (a) the result of QCD fit 1 (in which only quarks contribute to the pomeron structure at = 3 GeV^) 
is superimposed. In (b) the preferred QCD fit 3 (in which both quarks and gluons contribute) is shown. 
The results of the QCD fits are only shown in the kinematic range used in the fit (corresponding to data 
for which > 2GeV). In both figures the solid line represents the value of F^^^\xjp, /3,Q'^), whilst 
the dotted line shows the contribution from the pomeron term only. 
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Figure 8: The sum of the light quark distributions and the gluon distribution for fits 2 and 3, shown 
at (a) Q2 = 4.5 GeV^, (b) = i2GeV2 and (c) ^ 75GcV2. In fit 2 the gluon distribution is 
parameteriscd at the starting scale Qq = 3 GcV^ with only the first term in the polynomial expansion, 
while in fit 3 the first 3 terms are included. In both fits 3 polynomials are used to parameterise the 
light quark distribution. The figures are normalised such that they represent the parton distributions 
multiplied by the flux factor at x^p = 0.003. 
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Table 3: The quantity x^pF^^^^ as a function of /? and together with the statistical and systematic 
errors. Not included in the errors are separate overall normalisation uncertainties of ±6.0% for the data 
with <9 GeV^ and ±4.8% for those with > 9 GeV^. The data marked with an asterix correspond 
to the cross section, ^^"I^dQ^ ' averaged across the bin in /3. The /? boundaries for these data are 0.5 
and 0.8 for the bins at /3 = 0.65, and 0.8 and 1.0 for the bins at /3 = 0.9. 
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Table 3: The 
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x^Pp"^ (cont'd). 
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